#pragma once

#ifdef TXYZ

void print_propagator(colorcolorspinspin *S)
{
  printf("This is the propagator in space position \n");
  for (int t=0; t<T; t++)
        for (int x=0; x<L; x++)
          for (int y=0; y<L; y++)
             for (int z=0; z<L; z++)
             {
                 int dest=z+L*(y+L*(x+L*t));
                 for (int ic1=0; ic1<3; ic1++)
                   for (int ic2=0; ic2<3; ic2++){
                              for (int is1=0; is1<4; is1++)
                                   for (int is2=0; is2<4; is2++){
                                              printf("%8.16f\n",creal(S[dest][ic1][ic2][is1][is2]));
                                              printf("%8.16f\n",cimag(S[dest][ic1][ic2][is1][is2]));
                                                }
  
                        }
                }
}
#endif



#ifdef TYZX

void print_propagator(colorcolorspinspin *S)
{
  printf("This is the propagator in space position \n");
  for (int t=0; t<T; t++)
        for (int y=0; y<L; y++)
             for (int z=0; z<L; z++)
             for (int x=0; x<L; x++)
             {
                 int dest=z+L*(y+L*(x+L*t));
                 for (int ic1=0; ic1<3; ic1++)
                   for (int ic2=0; ic2<3; ic2++){
                              for (int is1=0; is1<4; is1++)
                                   for (int is2=0; is2<4; is2++){
                                              printf("%8.16f\n",creal(S[dest][ic1][ic2][is1][is2]));
                                              printf("%8.16f\n",cimag(S[dest][ic1][ic2][is1][is2]));
                                                }

                        }
                }
}
#endif
